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Abstract 



Based on the correspondence between Collins diffraction formula (optical Fresnel transform) 
and the transformation matrix element of a three-parameters two-mode squeezing operator 
in the entangled state representation (Opt. Lett. 31 (2006) 2622) we further explore the 
relationship between output field intensity determined by the Collins formula and the input 
field's probability distribution along an infinitely thin phase space strip both in spacial domain 
and frequency domain. The entangled Wigner function is introduced for recapitulating the 
result. 

OC1S codes: 070.2590, 270.6570 

In a preceding Letter [T] we have reported that the Collins diffraction formula in cylindrical 
coordinates is just the transformation matrix element of a three-parameter (fc and t are complex 
and satisfy the unimodularity condition kk* —tt* = X) two- mode squeezing operator [21 [3] 
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in the deduced entangled state representation (s, r'\ , 
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where ip s and 4> s denote the incoming and output light, respectively, [ai,a^] = 5i.j, 
k 



X - [A + D - i (B - C)} , t=^[A-D + i(B + C)} 
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we see that the relation kk* 



tt* = 1 becomes AD — BC = 1, J s is the sth Bessel function, and 
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here is the entangled states in two- mode Fock space [H[SJ[BJ[7] named after Einstein-Podolsky- 
Rosen (EPR)'s [5] concept of quantum entanglement, 



\V) = cx P {-i M 2 + v4 - V*4 + a\4} |00) 



(5) 



Thus (s, r'\ F^' k ^ is the quantum optics version of the Collins formula (generalized Hankel trans- 
formation). In [9] we have also found 
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Comparing with the integral kernel of usual Fresnel transform which describes how a general beam 
i\) (x') , propagating through an (ABCD) optical paraxial system, becomes output field <fi (x) [TUHTT] 
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where AD — BC = 1, 
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we see that K.^'^ (77', 77) can be considered as the integration kernel of 2-dimensional entangled 
optical Fresnel transform, 

*(t/) = fK^(r,',ri)9( V )d\ (9) 



in this sense F^'^ can be named entangled Fresnel operator (EFO), here <!> (77) = (ry| $) , ^ (77') = 
(77' I \P) , and we have used the completeness relation J 3-2 Jjy) = l. 

Clearly, if the [ABCD] system is changed to [13 (— B) (— C) A] system, then Eq. should read 



(10) 



where /C 



(D,-B,-C) . 



IS 



K {D,-B,-C) u 



1 



-2iB7T 



exp 



-2B 



D |?7| 2 — (7771'* + 77*77') + A |7/|' 



(11) 



On the other hand, signals or images in optical information theory may be described directly or 
indirectly by the Wigner distribution function (WDF) [12] . In one-dimensional (ID) case, the WDF 
of an optical signal field ip (x) is defined as 
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W^,{v, x) involves both spatial distribution information and space- frequency distribution information 
of the signal, v is named space frequency. Now, let us consider the entangled case. Like Eq. (|12p . 
it is natural to introduce the 2-D complex Wigner transform as 
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where a, 7, 77 are all complex variables. To see its physical meaning, using the integration formula 
of Dirac S— function, we perform the following integration, 
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which is just the probability distribution of the complex function 7/1(0"). Further, let the ordinary 
Fourier transforms of ip (a) be j (C) , 
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then substituting (| 1 5|) into (|T3|) leads to 
d 2 r/d 2 Cd 2 C . 
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It then follows from ([To]) that 



| dVW (a, 7 ) = J ^-j ( 7 - C) 3* (C + 7) j dW^* 



—3 (7 " 3* (t + l)6 (0 <J (O = - | j (7)| 2 : 

7T 7T 



(17) 



which is the probability distribution of the complex function j (7). Thus our definition in (flU)) leads 
to two marginal distributions in a and 7 phase space, respectively. Hence W (a, 7) is indeed the 
correct complex 2-D Wigner function ( Wigner transform) of complex function tp (a) or j (7) . If one 
wants to reconstruct the Wigner function by using various probability distributions, obviously the 
"position density" \(a\ ip}\ 2 and the space-frequency density K7I ip}\ 2 are not enough, so we extend 
S (77) S (77*) = 8 (771) 8 (772) to 8 (771 — Dai — Bj 2 ) 8 (772 — Da 2 + Bji) and generalize (fT4|) to, 
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R2 (771, 772) is also a probability distribution along an infinitely thin phase space strip denoted by the 
real parameters B,D, which is a generalized entangled Radon transform [13l [14] of the two- mode 
Wigner function (in the entangled form) [TSl US] , 

Then an interesting question naturally arises: what is the relation between the generalized Fresnel 
transform and the WDF in entangled state representation? 

We begin with rewriting the 2-D WF (JT3J) as 
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Substituting (fTO)) into |T5J) we rewrite the Radon transform of W (a, 7) as (d 2 cr = deride, d 2 7 
d7id7 2 ) 
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On the other hand, when the beam (77) propagates through the [D (-B) (— C) A] optical system, 
according to the Fresnel integration (fT0]) - (fTTj) . we have 
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which is the same as i?2 fai, 772) m (|2T)|) . So combining (|2H| . pT7)) - (fTT|) and (T2T|) we reach the conclusion 
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where AD — BC = 1. The physical meaning of Eq. ([22]) is: when an input field propagates through 
an optical \D(—B) (— C) A] system, the energy density of the output field is equal to the Radon 
transform of the two-mode entangled Wigner function of the input field. So far as our knowledge 
is concerned, this conclusion seems new. Eq. (|22[) is the relationship between the input amplitude 
and output one in spatial-domain. Next we turn to the frequency domain. 

If taking the matrix element of F^ ,fc ' in the |£) representation which is conjugate to I77), where 
the overlap (rj\ £) is (r/\ £) = | exp[(£?7* — £*7i)/2], we obtain the 2-dimensional GFT in its 'frequency 
domain', i.e., 
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where the superscript N of /Cj means that this transform kernel corresponds to the parameter 
matrix N = [D,—C,—B,A]. Thus if the [D, — C, —B, A] system is changed to N = [A, C, B, D] 
system, the GFT in its 'frequency domain' is given by 
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On the other hand, in similar to (|18p . we consider the integration transform, 

#2(6,6)=^ J 6^ 1 -Aa 1 -C l2 )S(^-Aa 2 +Cj 1 )W(a, 7 )d 2 ad 2 1 , (27) 

i?2 (6, 6) is also a probability distribution along an infinitely thin phase space strip denoted by the 
real parameters A, C. Substituting (fT§)) into (|2"T|) yields 
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This is the relationship between the output amplitude and input one's entangled Wigner function 
in 'frequency domain'. 

In sum, based on the correspondence between Collins diffraction formula (optical Fresnel trans- 
form) and the transformation matrix clement of a three-parameters two-mode squeezing operator in 
the entangled state representation, we have explored the relationship between output field intensity 
determined by the Collins formula and the input field's probability distribution along an infinitely 
thin phase space strip. The entangled Wigner function is introduced for recapitulating the result. 
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